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Abstract 

A discussion of how to calculate asymptotic expansions for polyhomogeneous spacetimes 
using the Newman-Penrose formalism is made. The existence of logarithmic Newman- 
Penrose constants for a general polyhomogeneous spacetime (i.e. a polyhomogeneous space- 
time such that "Jo = 0(r -3 ln^ 3 )) is addressed. It is found that these constants exist for the 
generic case. 



1 Introduction 

The study of the gravitational field of isolated sources close to null infinity J? has been done 
by solving the asymptotic characteristic initial value problem using formal series expansions in 
powers of 1/r, where the coordinate r is either the affine parameter of a future oriented null 
geodesic or a cleverly constructed luminosity parameter. However, as was noted by some of the 
first researchers in the area JTlJ , combinations of powers of In r and 1 jr tend to arise naturally, 
even if we assume initially expansions in 1/r. For long it was thought that these logarithmic 
terms were in some way connected to the presence of incoming radiation. In order to avoid the 
appearance of these logarithmic terms ad hoc, some conditions were imposed on the Ansatz for 
the asymptotic expansions. One of these conditions, used in the work by Bondi et al. §), @ 
and Sachs [^3| (Bondi-Sachs (BS) framework) precluded the appearance of 1/r 2 terms in some 
of the metric functions. This condition has been known for a long time as the outgoing radiation 
condition for the gravitational field (ORC), because it bears some resemblance to Sommerfeld's 
outgoing condition for the wave equation in flat space. It was later realized that the ORC does 
not preclude the existence of incoming radiation. Moreover, the presence of incoming radiation in 
itself is not a problem. There are some examples of spacetimes that contain a mixture of incoming 
and outgoing radiation and that have a somooth ' . In any case, what has to be avoided (and 
is implicit in the definition of an isolated body) is the presence of radiation travelling "infinitely 
long" distances (i.e. coming from past null infinity J^~). But this is achieved by imposing some 
boundary conditions directly on frL4| , 

Newman and Penrose, working with their spin coefficients formalism (NP formalism), used 
a different condition. Namely, they assumed the component of the Weyl tensor to behave 
asymptotically as 0(r~ 5 ). This condition has its roots in the conformal compactification tech- 
niques introduced by Penrose. If one assumes — 0(r -5 ) then one readily deduces the Peeling 
theorem, and all conformally rescaled quantities (the spin coefficients, the \I/'s and the tetrad 
functions) are well defined at J? . But if a slower decay for \& is assumed — say 0(r -4 ) — then 
logarithmic terms will appear again in the asymptotic expansions. One also finds that some 
conformally rescaled quantities diverge at ^ . Hence one ends up with a null infinity that is 
non-smooth. 
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Not surprisingly, the ORC and Penrose's condition are found to be related. If the ORC does 
not holds (i.e. if we retain the "nasty" 1/r 2 terms in our initial expansions) then H>q = 0(r~ 4 ) 

Hf. 

More recently, the work on the non-linear stability of Minkowski spacetime by Christodoulou 
and Klainerman has shown that generic spacetimes may exhibit a slower fall off (in some cases 
as low as 0(r -7 / 2 ) on ^o) and a non-smooth J? . These results together with the fact that 
logarithmic terms have been observed by other authors in the context of cylindrically symmetric 
spacetimes Q and infinite rotating discs [gj suggest that asymptotic expansions just in terms of 
powers of 1/r only are too restrictive for many physically interesting applications. 

From the previous discussion, one can see that a natural way to generalize (and hopefully gain 
more insight and better understanding!) the work that so far has been done on the asymptotic 
behaviour of the gravitational field is first, to assume a slower decay for the leading terms of the 
components of the Weyl tensor, and second, to allow the full appearance of logarithmic terms 
in the asymptotic expansions^. Spacetimes that can be expanded asymptotically in terms of a 
combination of powers of 1/r and In r are known as polyhomogeneous spacetimes. 

An analysis of the solutions of the constraint equations by Andersson & Chrusciel ||, Q 
seems also to indicate that an adequate framework to study the asymptotic gravitational field is 
that of the polyhomogeneous spacetimes. Chrusciel, MacCallum and Singleton have shown 
that the hypothesis of polyhomogeneity is formally consistent with the Einstein field equations. 
Many results that hold in the non-polyhomogeneous setting can be recovered. The BMS group is 
still the asymptotic symmetry group, and the uniqueness results for the Bondi mass of Chrusciel, 
Jezierski and MacCallum hold However, the well posedness and existence/uniqueness of the 
asymptotic characteristic initial value problem is still an open question. Other classical results 
like the Peeling theorem no longer hold, but an adequate generalization can be found in order 
to encompass the "new" logarithmic terms 

Among this last group of results one has that of the Newman-Penrose constants, a set of 10 
absolutely conserved quantities that was found to exist for spacetimes satisfying the = 0(r~ 5 ) 
condition. As shown in reference Jl(| in a particular example, the Newman-Penrose constants 
are no longer constants for polyhomogeneous spacetimes. However, it was possible to show that 
for a restricted class of polyhomogeneous spacetimes (those spacetimes with a finite conformally 
rescaled shear a at J^) it is possible to construct another set of 10 quantities (that we will call 
logarithmic Newman-Penrose constants) with a structure similar to that of the NP constants and 
that are conserved [ pc| . 

The objective of this article is to extend the work of references Jl(| and |2(| on these logarithmic 
constants, i.e. to find a family of functionals of that for an arbitrary polyhomogeneous 
spacetime yields conserved quantities, and for a non-polyhomogeneous spacetime reduces to the 
NP constants. In order to do so, it will be necessary to study in some detail the polyhomogeneous 
expansions of a spacetime. Some acquaintance with the NP spin coefficients formalism will be 
assumed in order to follow the calculations. 

Much has been discussed about the interpretation and physical meaning of the NP constants: 
they have been regarded by some authors as the gravitational analogue of the electric charge, as 
they retain their values even after a process of radiation ]16[ |, |l8[ . However, the works trying 
to find an interpretation for them have been so far inconclusive, and some authors have even 
suggested on some grounds that the constants may be devoid of physical meaning [Q , |2^] , |l^] . 
The new set of logarithmic constants adds to this old debate. Once their existence has been 
shown, the next natural step is to try to attach to them a physical meaning (if any!). In order 
to do a connection with previous work, stationary polyhomogeneous spacetimes are studied, and 
their logarithmic NP constants calculated. The result shows that for stationary spacetimes the 
logarithmic constants have a "logarithmic quadrupole" structure. 

This article is organized as follows: in section 2 we give a brief description of the notations 
and conventions that are used in the present aerticle. In section 3 there is a discussion on how 
to calculate asymptotic expansion for polyhomogeneous spacetimes using the Newman-Penrose 
formalism. A description of the tetrad choice is also given. In section 4 we address the problem 

1 A precision should be done here. We will exclude from our study those spacetimes that contain infinite series 
in powers of ln(r). 
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of how to generalize the Newman-Penrose constants to polyhomogeneous spacetimes (logarithmic 
NP constants). A detailed study of the evolution equation for the component of the Weyl 
tensor is done. Using this analysis, the main result of this article is proven: the existence of the so 
called logarithmic Newman-Penrose constants for all polyhomogeneous spacetimes. Finally, there 
is one appendix. In it one finds 2 tables that summarize the behaviour of generic polyhomogeneous 
spacetimes. The results listed in these tables are used extensively in most of the arguments in 
this paper. 

2 Notation & conventions. 

We will be working with spacetimes that formally can be expanded in terms of powers of 1 /r and 
lnr. For example 

N k 

o = ^Ta k r- k = ^^o kl r- k \n> r, (1) 

k k j 

where Uk denotes the polynomial in z = In?' with coefficients depending on u, 8 and <p associated 
with r~ k . A second subscript will mean the power of z to which this function is attached. So, aij 
is the coefficient that comes with r~ l In- 7 r. For the components of the Weyl tensor superscripts 
rather than subscripts will be used. Let the symbol # denote the degree of a polynomial in z. 
We will stick to the convention that the zero polynomial has degree — oo. Hence 

#c(u, 0,^0 = 0, (2) 

#0 = -oo. (3) 

A number in square brackets to the right of a polynomial will also indicate its degree. For instance: 
a 2 [N 3 + 1] means that #er 2 = N 3 + 1. 

We will refer to the NP "field equations", the Bianchi identities and the frame equations in 
the way they are listed in Stewart's book [Q. In many ocasions one will only be interested in the 
terms that go with a determinate power of 1/r. This will be indicated by attaching a subindcx 
to the name of the equation. A second subindex will refer to a particular power of lnr. So, for 
instance (Bb^jVg will mean that we are looking at the relationship that goes with the r 7 ln^ 6 
term of the Bianchi identity Bb. 

When calculating the expansions of the diverse quantities, one usually needs to evaluate the 
derivatives and products of polyhomogeneous functions. The following formulae are most useful 
in these circumstances: 

hg = J2J2 h k( z )9i-k(z)r-\ (4) 

i=2 k=l 

= E - (i ~ 1) hi-x{zj) r-\ (5) 

i=1 

where h and g are two arbitrary polyhomogeneous functions, the coefficients g k are polyno- 
mials in z — lnr, and ' denotes differentiation with respect to z. 

Some use of the spin weighted spherical harmonics ( s Yi m ) and the 3 and 9 (eth and ethbar) 
operators will be made. In this case, we will also stick to Stewart's convention that is the same 
as the one in Penrose and Rindler Ell . Their connection with the NP directional derivatives is 



= 6 + 8(3-0), (6) 
= 6-s(a-]3). (7) 
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3 Solving the NP equations using polyhomogeneous expan- 
sions. 



The first step into solving the NP field equations is to construct a coordinate system and a null 
tetrad adapted to the problem. Suppose there is a one parameter family of null hypersurfaces. 
Denote this family by Af u . We can use the parameter u to define a scalar field. Of course, this 
paramctrization is completely arbitrary, and one could have used another parameter u as long as 
some smoothness requirements are satisfied. On each of the geodesic generators j u of the null 
hypersurfaces one can choose an affine parameter r. The null hypersurfaces M u intersect in 
cuts S u . On any one of the cuts, say So we can choose arbitrary coordinates x a , a = 2, 3. We 
can propagate the coordinates x a on J 1 by demanding x a — const, on each generator of J? . In 
this way we have constructed a coordinate system (u,r, x a ) in the neighborhood of J? . There is 
some freedom left in this construction: 

1. we can always choose a different parameter u for the family of null hypersurfaces; 

2. we can choose a different selection of coordinates x a on Sq,x" — x Q (a; a ); 

3. we can reset the origin and scaling of the affine parameter r, r = a(u, x a )r + b(u, x a ). 

Once we have set our coordinate system, we can construct a null tetrad. An obvious choice 
for the first element of the tetrad is the gradient of the null hypersurfaces 



1 = du. 



(8) 



The vector field is tangent to the geodesic generators 7„ of Af u . We can use the freedom of the 
scaling of the affine parameter of these generators to set 



dr 



(9) 



The freedom of the choice of the origin of the affine parameter will be used so that p will not 
contain any 1/r 2 terms (see equation (p^)). Let us denote the 2-surfaces u — const., r — const. 
by S Utr . It can be seen that is future pointing and orthogonal to S u _ r . There is only one 
other null direction with the same properties. The second vector of our tetrad n M is chosen to 
be parallel to this direction. The two other vectors of the tetrad, m M and in' 1 , are constructed so 
that they span the tangent space to <S UiI ., T(S Ujr ). So from this construction one deduces 



d d d 
ou or ox a 



(10) 



so that the contravariant metric tensor is 



? dx a ' 



(11) 



{lis 

9np 



(12) 



/ 1 \ 

1 2Q C 9 C v 

C 6 -2i e f 

\ C v - £V -2^f j 

Using a spin-boost it is possible to set e = e everywhere. From the commutator equations we 
deduce 



K = e = 0, 



(13) 



a + /3 



(14) 



and p and p are real. 
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The asymptotic characteristic initial value problem is usually set by supplying "Jo on an initial 
null hypersurface Af , a on J^+, and *i, Re* 2 and f on 2 = J+ n7V Q. 

The most general polyhomogeneous form for that is physically reasonable to start with is 

*o = X>o, (15) 

k=l 

where the coefficients are polynomials in lnr, as mentioned before. With this initial data, the 
first two equations in the NP hierarchy to be solved are 

Dp = p 2 + aa, (16) 
Da = 2pa + f . (17) 

Using (|) and ^ into @ and @ one obtains the following recurrence relations 

i-l 

Pi-l - (* _ 1) Pi-l = X! (PkPi-k + CfcOVfc) , (18) 

fc=l 

i-l 

- (* - 1) o-i-x = 2 ^ p fe( 7 4 _ fc + (19) 
fc=l 

Setting i = 1 we find that 

vfo = 0. (20) 
For i = 2, the recurrence relations reduce to 

Pi ~Pi=Pi + > ( 21 ) 
- a x = 2 Pl a x + *g. (22) 

From equation (|2l| ) we see that pi has to be a polynomial of degree zero. Couch & Torrence Q 
have shown that in order to have a spacetime with a Weyl tensor that vanishes at infinity then 
we require \l/o = 0(r -2-£1 ) where e± > 0. Hence, in order to have a physically realistic situation 
we set 

*o = 0, (23) 
so that a\ is also of degree zero. The remaining system then reduces to 

-px =Pi + axai, (24) 
-ai = 2 Pl ax. (25) 

which has 3 solutions: p\ = 0, a\ = (asymptotically plane) which yields degenerate spacetimcs 
for which = m M = 0; p\ = —\ and \ai\ 2 = i (asymptotically cylindrical) which also produces 
degenerate spacetimes for which — m M ; and p\ = — 1, a\ = which gives rise to spacetimes 
that are asymptotically minkowskian. In what follows we will stick to this last solution. 
For i = 3 one finds that 

P2 = 0, (26) 
a' 2 = *g. (27) 

Therefore p 2 is a polynomial of degree in lnr. Using the remaining freedom in the definition 
of the radial coordinate (recall that r is an affine parameter) it is possible to redefine r so that 
Pi = (see for example |ig[|). From equation (p7|) one reads 

^ = (j + l)a 2j+1 , (28) 

j = . . .N3, whence #02 = N3 + I. So, the only part of 02 that is not determined by the data on 
the initial hypersurface TV is 02.01 i n complete agreement with the results on the well posedness 
of the asymptotic characteristic initial value problem. We can think of 'J'q'- 7 j — 0...N — 1 as 
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a contribution to the shear from the incoming radiation using Szekeres interpretation of the 
components of the Weyl tensor (\l/o can be regarded as an incoming transverse wave) [B5| . 
For i = 4 one obtains, 

P3 - P3 = (29) 
4 - a 3 = *S. (30) 

From these two equations one can readily see that = 2 A3 + 2 and #173 = N4. Recurrence 
relations to calculate the coefficients of the two polynomials can be easily obtained. 

The process described here can be performed up to any desired order in 1/r. A similar analysis 
can be done with the remaining field equations, Bianchi identities and commutator relations. The 
order of solving the equations is given in reference [ p4| . A summary of this analysis is presented 
in the form of tables in Appendix A. 



4 Logarithmic NP constants for a generic polyhomogeneous 
spacetime. 

We are interested in the construction of conserved quantities that can be expressed as integrals 
over a cut of J? , S 2 . One expects these constants to be formed out of quantities that are initial 
data on the initial null hypersurface Afo- As "J^' , Re^' , and 02,0 (the data that has to be 
prescribed on J? + n Afo) are used to construct the Bondi mass and the angular momentum of 
the spacetime — quantities that are not conserved — , we do not expect to be able to construct 
conserved quantities out of them. 

Newman & Penrose jl(|, |l8| have found that for non-poly homogeneous spacetimes, 10 con- 
served quantities can be constructed out of ^q. As it has been mentioned before, these quantities 
are not conserved for polyhomogeneous spacetimes. Nevertheless, in the case when the leading 
term of the shear (02) contains no Inr terms (which corresponds to the situation = —00) 
it was possible to single out some other quantities with a similar structure (logarithmic NP 
constants) that are conserved. We want now to show that even for the most general class of 
polyhomogeneous spacetimes it is possible to construct such constants. 

The key idea in the construction of these constants is to deduce from the Bianchi identities 
an equation whose left hand side is a derivative with respect to the retarded time u of one of 
the coefficients of ^0, and whose right hand side is an S derivative of some combination of the 
coefficients in ^0, ^2; that is, equations of the form 

^• M =s(n (3i) 

where ' denotes differentiation with respect to u, and the indices i and M are to be determined. 
Note that this equation has the same structure as the continuity equations in mechanics of 
continuous media: a derivative with respect to time plus the divergence of a flux, F Equation 
( pl| ) has an overall spin weight 2, and the function F spin weight 3. 

We multiply equation ( |3l| ) by {vYim), a spin 2 weighted spherical harmonic, in order to obtain 
an equation of spin weight 0. Now, we integrate over the unit sphere (a cut of so that 



d_ 

Oil 



^6 M (2Y lm )dS) = I B(F)( 2 Y lm )dS. (32) 



Using the identity E( 



S 2 



.Yt, m )ti S+1 CdS = 0, (33) 



l S 2 

where C is of spin weight I + 1 (in our case C = F has spin weight 3, hence 1 = 2), we end up with 



2 The operator 93 applied to a quantity of spin weight zero is equal to the angular part of the laplacian, V 2 . 
Hence, we can regard 99 as the laplacian for spin weighted quantities, and 9 or 9 as divergences ("square root of 
the laplacian"). 
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0_ 

du \J S 2 



n M (2Y 2m )dS ) = 0. (34) 



Hence the 5 complex quantities 



Ql ,M = / ^M {2 Y 2m )&S, (35) 
Js 2 

m = —2, —1 ... 2 are conserved as long as equation (|3l| ) holds. 

The results in |^6| show that one should expect to find the desired conservation laws in the 
coefficients of the highest power of lnr in the 1/r 6 terms of the evolution equation for <3>o- An 
algebraic explanation of this fact can be found by looking at the identity 

(A(u, 0, tp)r- 1 ln J r) = {-i)A{u, 6, tp)^' 1 W r + jA{u, 9, p)r _i_1 W 1 r, (36) 

Hence, differentiation with respect to r adds extra terms to the lower powers of lnr, but not to the 
highest. Therefore, the coefficients of the highest power of lnr (for a given power of 1/r) mimic 
the behaviour of the coefficients in a non-polyhomogeneous spacetime. The evolution equations 
in both cases will have the same algebraic form. 

The situation discussed in this article is more elaborate than the one described in [p6| as 
one has to deal with many new terms coming from 'J'q. It will be necessary to keep track of 
the behaviour of the degrees (Aj) of the polynomials as the system evolves. Therefore our 
discussion will start at the lowest order (1/r 3 ) instead of just going directly to the 1/r 6 order 
terms. 

As discussed previously, the most general form for ty is the following: 

* = *g[JV 3 ]r- 3 + *Mr- 4 + * 5 [N 5 ]r- 5 + ^[A 6 ]f- 6 + • ■ • , (37) 
and its evolution equation is the Bianchi identity (Bb), 

A* - = (47 - /x)* - 2(2r + + 3aV 2 - (38) 

One expects the values of Ni to change as the system evolves. So, we can think of these numbers 
as discrete functions of the retarded time u (Ni — A* (it)). A tilde ~ over one of these numbers 
will mean that we are referring to its value at the initial null hypersurface A/o- 
The equation (Bb)3 readily yields 

* g = 0, (39) 

whence the polynomial VP 3 , does not depend on u, and its degree (#'I ,3 ) = -^3 = -^3) ^ s preserved 
during evolution. 



From (Bb)4 we obtain 



2 

an equation of the form 



*o - ^o'[^3 - 1] + *g[JVs] - 3*? [A 3 ] = 0, (40) 



*o = W], (41) 



P[Ns] a polynomial of degree A3. If at the initial null hypersurface A/q we have N4 > A3 then 
we find that 

= = . . . = ^N a +1 = Qj (42) 

hence ^q' N,l , . . . , ^q' 3+1 will be constants of motion, and the number A4 a constant. The first 
evolution equation where the interaction between r -3 and r -4 terms shows is 

Using, (Ba) 4 

* 3 ' [ A 3 - 1] + * 3 [A 3 ] = S* 3 [N 3 ] , (44) 
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and the commutator for the 9 derivatives we obtain 

^ Na + 3^ Ns = m^ 3 - Ns , (45) 

so that we cannot derive a conservation law from this equation due to the presence of the extra 
3*q term. If 7V 3 > N A then the only thing we can affirm is that iV 3 > 7V 4 (from (|l|)). 
The (Bb)5 equation is much more involved. 

*o[^s] - - 4^) + Qi[JV 3 + 1](M# - 3*3) [N 3 ] + C«{N 3 + l]d a *l[N 3 ] 

-S^i[Y] + 5 2 [N 3 + l]^l[N 3 ] = 
4 72 [7V 3 + l]*g[JV 3 ] + i*J[JV4] - fi 2 [N 3 + l]*g[JV 3 ] - 4r 2 [7V 3 + 1]*?[JV 3 ] 

-2/3 x *?[Y] - 2f3 2 [N 3 + 1]*?[JV 3 ] + 3a 2 [iV 3 + 1]*|[JV 3 + 1], (46) 

an equation of the form 

*o = P[Y] (47) 

where Y = max{iV 4 + 1, 2A" 3 + 2}. Note that if iV 4 + 1 > 2N 3 + 2 then N 4 > N 3 , and again 
due to (Bb) 4 , 7V 4 will be a constant. Hence in this case Y is a constant. If in the other hand, 
N 4 + 1 < 2N 3 + 2 we end up with two cases (N 3 < N 4 or N 3 > N 4 ). If N 3 < N 4 then N 4 is a 
constant, and_then the inequality 2iV 3 + 2 > N 4 + 1 holds for later times because it held at A/o- 
And if A" 3 > N 4 then N 3 > N 4 , and hence again 2N 3 + 2 > N 4 + 1. The conclusion is that Y is 
a constant again. 

Looking at equation (Eft) we find that if 7V 5 > max{A^ 4 + 1,2AT 3 + 2} we obtain a set of 
constants of motion (*q' , . . . , ^ Y+1 ). If A5 < Y = m&x{N 4 + 1, 2N 3 + 2} at Af then we 
will not have constants of motion. The first equation showing an interaction between terms of 
different orders in r is equation (Bb)5.y: 

^ Y = Ml> Y + {* 2 , N3+1 *l N > +1 }, (48) 

where the term in curly brackets {} is present if Y = 2N 3 + 2. In order to obtain an equation in 
terms of "initial data" quantities only we require the equations (Ba)s and (Bc) 4 : 

*?[Ni] = m 4 [N 4 ] + {f - Aa 2 )[N 3 + 1}V 3 [N 3 ] + ir 2 [N 3 + 1}V 3 [N 3 }. (49) 

* 3 '[N 3 ] -mi[N 3 ] = -\t[0]mN 3 ] (50) 

Using equations (^) and (|5(]) it is not hard to show that equation (Q) does not have the required 
form. 

Following jl8| and p(| we expect to find conserved quantities of the desired form (if any) in 
the (Bb)g terms. The equation in this case is extremely involved, so it is better to look at it in 
parts: 
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4(*o'-5*o) 



-Qi[JV 3 + l] (*o'-4*o) [Na] 



(51) 





+Q 2 [Y] (^'-3^) [jv 3 ], 


(52) 




= C?[N 3 + l}d a *t[N 4 } + C%[Y]d a V 3 [N 3 }, 


(53) 




= S 1 ^ 5 1 [N 5 ]+S 2 [N 3 + l]^i[Y] 






+S 3 [max{N 4 ,2N 3 + 2}}^ 3 [N 3 }, 


(54) 




= 72 [ N 3 + 1] ^ [jv 4 ] + 73 [ Y ] *s [JV 3 ] , 


(55) 




= -5*0^5] + m[N 3 + 1]*J[JV4] + » 3 [Y]V 3 [N 3 ], 


(56) 


[^i] 6 


= T 2 [N 3 + l}*i[Y]+T 3 [Y]H> 3 [N 3 ], 


(57) 




= /3^l[N 5 ] + /3 2 [N 3 + l}^f[Y} + I3 3 [Y]^ 3 [N 3 ], 


(58) 




= a 2 [N 3 + + a 3 [N 4 }^ 3 [N 3 + 1]. 


(59) 
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Hence, (Bb)6 is an equation of the form 

*o = P[X] (60) 

where X = max{iV 5 , 3iV 3 + 3, N 3 + N 4 + 2} = max{jV 5 , (N 3 + 1) + Y}. Now, if initially N 5 > 
(N 3 + 1) + Y then N$ > Y, and using (Bb)s we see that N$ is a constant. On the other hand, 
if N5 < (N 3 + 1) + Y then we end up with two possibilities: N$ > Y which yields N$ constant; 
and N5 < Y from which one deduces N5 < Y. This arguments show the remarkable fact that the 
number 

X = m&x{N 5 ,3N 3 + 3,N 3 + N 4 + 2} (61) 

is a constant independent of the initial values of N 3 , N4 and N5. 

As in previous orders, we will have a set of constants of motion if N& > X = ma,x{N^, 3-/V3 + 
3, N 3 + N 4 + 2}: ^ Nli , . . . , The first non-trivial evolution equation will be 



+ {3*2,n 3 +i*% Y } . (62) 

The terms in curly brackets { } will not be present if X = N5, while the terms in angle brackets 
( ) will not appear if X = max{3iV 3 + 3, JV 3 + i\T 4 + 2}. 

In order to obtain an equation that contains only "initial data" quantities we will need to use 
the equations (Ba)6 and (Bc)s: 



*f [N 5 - 1] - **[N 6 ] - m 5 [N 5 ] - f[N 3 + 1]*§[JV 3 ] 
-f[max{N 4 2N 3 + 2}}^l[N 3 ] = 
MN3 + iMiNi] + n 3 [Y]^ 3 [N 3 ] - 4a 2 [N 3 + 1]*S[JV 4 ] 

-4a 3 [Y]* 3 [N 3 ] + 4p 3 [2N 3 + 2}* 3 [N 3 ], (63) 



-2 r 



^[Y - 1] - - m\[Y] - 6 [N 3 + IM[N 3 ] = 

-Ai*S[JV 4 ] - A 2 [iV 3 + l]*g[JV 3 ] + 2n 2 [N 3 + 1]*?[JV 3 ] - 2a 2 [N 3 + 1]#?[JV 3 



From equation (p3|) we obtain (when X = N 5 ) 



and from 



Substitution of these results into equation (62) yields 

33*^) - (2*^) - {(<5 2 ^ +1 - 2/3 2)JV3+1 )^' r } 
{4^3+1^} - {3<7 2lJ v 3+ i3^ y } • 



4# x 



(64) 

(65) 
(66) 



(67) 



Using the techniques described in the previous section is not difficult (but very lenghty!) to obtain 
the values of the coefficients 5 2 ' N3+1 , /?2,jv 3 +i, t 2,n 3 +i, and C2,jv 3 +i in terms of ^q: 



r2,jv 3 +i 



(68) 
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= tt^-^ 3 , (69) 



^3 + 1 

T2,n 3+ i = A^Vl^o'^ 3 . (70) 
= ^Vl^o^- (71) 



With them and the commutator relation 



(95 - 99)4-^ = 2<$l' x , (72) 



one obtains an equation with the same structure as (|3l"l): 

= -9 (jw**) + {^r^r^r}) (73) 

At this point, the relevance of the constancy of X becomes clear. If it were not the case, then 
the continuity equation ((73]) could change as the system evolves, and we would have different 
constants at different null hypersurfaces! 

Therefore we have proven the following proposition: 

Proposition 1 The quantities (5 complex in general) 

Qk=[ * 6 ' X ( 2 Y 2k )dS (74) 
Js 2 

with X — max{iV5, 3-/V3 + 3, -/V3 + N4 + 2} are constants of motion for any poly homogeneous 
spacetime. 

Observe that if N 3 = N4 = — 00 and N 5 = (that encompasses the non-polyhomogeneous 
case studied by Newman and Penrose) then X = 0, and the logarithmic constants reduce to 
the NP constants. While if N 3 — —00 then X = N 5 and we recover the logarithmic constants 
of fjFjf . Hence we have succeded in our primary objective of finding a family of functionals 
of the quantity 'J'o that are absolutely conserved, and that reduce to the NP constants in the 
non-polyhomogeneous framework. 



5 Conclusions 

It has been shown that there is a functional of that for a polyhomogeneous spacetime yields 
10 absolutely conserved quantities, and for a non-polyhomogeneous spacetime reduce to the NP 
constants. 

There are still some remaining open questions. One of them is related with the well known 
fact that for linear fields one has an infinite hierarchy of NP constants. The general feeling in the 
case of the NP and logarithmic NP constants for the gravitational field is that the only quantities 
that one can form that have the desired properties are precisely the ones already known. However 
a proof is still not available. 

The other main question is that of the physical interpretation of the quantities. A step forward 
in this line would be the evaluation of the constants for some non-stationary spacetime. These 
topics will be the subject of future work. 
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A Leading terms for the spin coefficients 



The objective of the present appendix is to present in a concise way the results of the calculations 
described in section 3. The tables give the degrees of the different polynomials that appear 
in the polyhomogeneous expansions. These tables are indispensable if one wishes to follow the 
calculations in full detail. As mentioned before, a — oo means that the given polynomial is zero. 

A.l General polyhomogeneous spacetimes (N 3 ^ — oo). 





l 


r- 1 


T-' 1 


r -3 


r- 4 




— oo 


— oo 


— OO 


N 3 


N 4 


*1 


— oo 


— oo 


— oo 


N 3 


max{7V 4 + 1,2-ZV 3 + 2} 




—oo 


— oo 


— oo 


N 3 + l 


max{7V 4 + 1,2AT 3 + 2} 


*3 


—oo 


— oo 









*4 


—oo 












P 


—oo 





— oo 


2N 3 + 2 


N 3 + N 4 + 1 


a 


—oo 


— oo 


N 3 + l 


N 4 




a 


—oo 





N 3 + l 


max{7V 4 + l,2iV 3 + 2} 




P 


—oo 





N 3 + l 


max{N 4 + 1,27V 3 + 2} 




T 


—oo 


— oo 


N 3 + l 


max{7V 4 + 1,2JV 3 + 2} 




7T 


—oo 


— oo 


N 3 + l 


max{7V 4 + l,2JV 3 + 2} 




7 


—oo 


— oo 


N 3 + l 


max{7V 4 + l,2N 3 + 2} 




A 


—oo 





N 3 + l 


max{iV 4 + 1,2JV 3 + 2} 




P- 


—oo 





N 3 + l 


max{JV4 + l,2JV 3 + 2} 




V 


—oo 


— oo 


N 3 + l 


max{iV 4 + l,2iV 3 + 2} 






—oo 





N 3 + l 


max{Ar 4 , 2JV 3 + 2} 




c a 


—oo 


— oo 


N 3 + l 


max{7V 4 + 1,2JV 3 + 2} 




Q 





N 3 + l 


m&x{N 4 + 1,2N 3 + 2} 







A. 2 Polyhomogeneous space-times with a finite shear at J . 





l 


r -l 


r~' 2 


r -3 


r -4 


r -b 


r -tt 




— oo 


— oo 


— oo 


— OO 


N 4 


N 5 




*1 


— oo 


— oo 


— oo 


— oo 


N 4 + l 


N 5 






— oo 


— oo 


— oo 





N 4 + l 


N 5 




*3 


— oo 


— oo 














\I>4 


— oo 
















P 


— oo 





— oo 





N 4 






a 


— oo 


— oo 





iV 4 


N 5 






a 


— oo 








N 4 + l 








P 


— oo 








N 4 + l 








T 


— oo 


— oo 





N 4 + l 








7T 


— oo 


— oo 





N 4 + l 








7 


— oo 


— oo 





N 4 + l 








A 


— oo 








A^4 + l 








M 


— oo 








A^ 4 + l 








V 


— oo 


— oo 





N 4 + l 








r 


— oo 








N 4 








c a 


— oo 


— oo 





N 4 + l 








Q 








N 4 + l 
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